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We present detailed numerical calculations of the light element abundances synthesized in a Uni- 
verse consisting of matter- and antimatter- domains, as predicted to arise in some electroweak 
baryogenesis scenarios. In our simulations all relevant physical effects, such as baryon-antibaryon 
annihilations, production of secondary particles during annihilations, baryon diffusion, and hydro- 
dynamic processes are coupled to the nuclear reaction network. We identify two dominant effects, 
according to the typical spatial dimensions of the domains. Small antimatter domains are dissi- 
pated via neutron diffusion prior to 4 He synthesis at T4 Hc « 80 keV, leading to a suppression of 
the primordial 4 He mass fraction. Larger domains are dissipated below Ti He via a combination of 
proton diffusion and hydrodynamic expansion. In this case the strongest effects on the elemental 
abundances are due to p 4 He annihilations, leading to an overproduction of 3 He relative to 2 H and 
to overproduction of 6 Li via non-thermal nuclear reactions. Both effects may result in light element 
abundances deviating substantially from the standard Big Bang Nucleosynthesis yields and from 
the observationally inferred values. This allows us to derive stringent constraints on the antimatter 
parameters. For some combinations of the parameters, one may obtain both, low 2 H and low 4 He, 
at a common value of the cosmic baryon density, a result seemingly favored by current observational 
data. 



26.35.+c,98.80Cq,25.43.+t 
I. INTRODUCTION 

Recently, there has been a revived interest in antimat- 
ter cosmologies, stimulated in part by the first flight of 
the Alpha Magnetic Spectrometer (AMS) and by the 
prospect of a long-term AMS mission on board of the 
International Space Station Alpha (ISSA). While these 
efforts concentrate on direct detection of antimiclci in 
the solar system, several limits on the existence of an- 
timatter domains have been placed in the past. It has 
been known for a long time, that the presence of signifi- 
cant amounts of antimatter within a distance of about 20 
Mpc from the solar system may be excluded on grounds 
of the non-observation of annihilation radiation Q . More 
recent studies of the diffuse 7-ray background claim to 
exclude antimatter regions in today's Universe within a 
distance of ~ 1000 Mpc g, a considerable fraction of the 
present horizon ~ 3000 Mpc. The existence of antimat- 
ter domains would also have impact on the Cosmic Mi- 
crowave Background Radiation (CMBR). Recent studies 
predict 'ribbon'- or 'scar'-like anisotropies in the CMBR, 
at the interfaces of matter- and antimatter domains, with 
a Sunyaev-Zel'dovich-type y distortion of the order of 
10~ 6 HQ- Such small distortions are beyond the detec- 
tion limits of current CMBR observations, and probably 
also beyond those of the upcoming MAP and PLANCK 
satellite missions. Given that only a small region of the 
parameter space for the sizes of antimatter domains re- 
mains, it seems very unlikely that we live in a Universe 
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containing any considerable amount of antimatter today. 
In particular, the possibility of a baryo-symmetric Uni- 
verse containing equal amounts of matter and antimatter 
is excluded unless the separation length of matter and 
antimatter is nearly as large as the current horizon. It 
is, however, not possible on grounds of the above results 
to exclude the existence of small and distant pockets of 
antimatter M. 

A complementary scenario with small scale domains 
of antimatter which have completely annihilated prior to 
recombination has, however, hardly been investigated in 
the past. Note that such a scenario necessarily involves 
an excess of matter over antimatter. While the very pre- 
cise observation of the CMBR allows us to place stringent 
limits on any non-thermal energy input into the CMBR 
during epochs with CMBR temperature 0.3 eV < T < 1 
keV, in particular, also on energy input due to annihila- 
tions H , even more stringent constraints may be derived 
from considerations of Big Bang Nucleosynthesis (BBN, 
hereafter) . Such a baryo-asymmetric Universe filled with 
a distribution of small-scale regions of matter or antimat- 
ter may, for example, arise during an epoch of baryoge- 
nesis at the electroweak scale. It has been shown within 
the minimal supersymmetric standard model, and under 
the assumption of explicit as well as spontaneous CP vi- 
olation, that during a first-order electroweak phase tran- 
sition the baryogenesis process may result in individual 
bubbles containing either net baryon number, or net anti- 
baryon number M. More recently, it has been argued 
that pre-existing stochastic (hyper)magnetic fields in the 
early Universe, in conjunction with an era of electroweak 
baryogenesis, may also cause the production of regions 
containing either matter or antimatter Though it 
seems questionable at present if electroweak baryogene- 



1 



sis has occurred at all, there are other imaginable scenar- 
ios which could lead to a small-scale matter- antimatter 
domain structure in the early Universe ||. 

This kind of initial conditions may have profound con- 
sequences on the abundances of the cosmological synthe- 
sized light elements. In the standard picture, synthesis of 
the light elements takes place between the cosmological 
epoch of weak freeze out at T w 1 MeV and T w 20 keV. 
The abundances of the light elements are highly sensi- 
tive to the cosmic conditions during that epoch. For 
example, the primordially synthesized 4 He mass fraction 
is sensitively dependent on the relative abundances of 
protons and neutrons at T4 Hc « 80 keV, when prac- 
tically all available neutrons are incorporated into 4 He 
nuclei. We have recently shown that annihilation of an- 
timatter domains during BBN (at temperatures above 
T4 Hc ) may significantly alter the neutron-to-proton ratio 
at T4 Hc pO| j. However, light element abundances are also 
sensitive to putative matter-antimatter annihilations af- 
ter the epoch of BBN. In this paper, we thus extend our 
analysis of annihilation of antimatter domains to a much 
wider temperature regime, from above the epoch of weak 
freeze out to the epoch of recombination. This allows us 
to constrain matter-antimatter domains within a much 
wider range of domain separations. The same scenario, 
annihilation of matter-antimatter domains during and af- 
ter BBN, has been very recently investigated in a Letter 
by Kurki-Suonio and Shivola. Jll]]. Though the main 
conclusions of our paper are not vastly different from 
those of Ref. jll] , we arrive at somewhat different results 
(factor ~ 3) for the synthesis of some of the elemental 
abundances. Furthermore, we use a different approach 
in comparing our theoretical results with observationally 
inferred values for the light element abundances. In par- 
ticular, we base our constraints on observationally deter- 
mined limits on the primordial 3 He/ 2 H ratio, rather than 
on the much less secure limit on primordial 3 He. We also 
consider the production of 6 Li which yields a tentatively 
much more stringent limit on the existence of antimatter 
domains. 

Prior studies of the influence of antimatter domains on 
the light element abundances have only been carried out 
in the context of a baryo-symmetric Universe E 12 1. Of 
course, such models have to assume that annihilation of 
all cosmic baryons may be avoided by an assumed "un- 
physical" , and unknown rapid separation mechanism of 
matter from antimatter. In essence, these works have 
shown that antimatter domains and successful BBN mu- 
tually exclude each other in baryo-symmetric cosmologies 
unless the separation between matter and antimatter do- 
mains is exceedingly large. In that case, however, BBN 
proceeds in a standard way, independently in matter and 
antimatter domains. 

There have been a number of studies concerning a ho- 
mogeneous injection of antimatter into the primordial 
plasma during, or after BBN ]T^-[T^|. Antibaryon pro- 
duction may result through the decay of relic, heavy par- 
ticles X, if hadronic decay channels are present, or the 



evaporation of primordial black holes. A possible candi- 
date for the A-particle is the gravitino, the superpartner 
of the graviton. It was realized early that antibaryons 
injected around weak freeze out would increase the syn- 
thesized 4 He abundance, due to preferential annihilation 
on protons 

0- This was used to derive a limit on the 
relative antibaryon abundance of ni/(n n + n p ) < 1/20. 
Injection of antibaryons after BBN would result in p 4 He 
annihilations The concomitant production of 2 H, 

3 H, and 3 He, either by direct production during 4 He an- 
nihilation, or by fusion processes of secondary neutrons, 
was found to give unacceptable large ( 2 H + 3 He)/H, un- 
less n^/nt, < 10~ 3 Q. Such arguments were subse- 
quently used to constrain the abundances of particular 
relic particles j^]. Motivated by the idea to reconcile a 
Universe dominated by baryonic dark matter with BBN, 
Yepes & Dommguez-Tenreiro investigated the effects of 
antibaryons injected during BBN Jf6| . An outstanding 
feature of such scenarios is a significant reduction of 7 Li 
production compared to standard BBN. Nevertheless, 
the claim that such scenarios may be compatible with 
a fractional contribution of baryons to the critical den- 
sity of Qb = 1 seems not viable due to the overproduction 
of 4 He and the 3 He/ 2 H ratio. 

In this paper we investigate BBN with matter- 
antimatter domains. Assumed initial conditions and def- 
initions which are used throughout the paper are intro- 
duced in Sec ||. In contrast to homogeneous antimatter 
injection, this topic requires an understanding of hydro- 
dynamic and diffusive processes which lead to the mixing 
of matter an d a ntimatter. These processes are summa- 
rized in Sec. 
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Nuclear annihilation reactions and the 
production of secondaries and their evolution are inves- 
tigated in Sec. IV. In Sec. [v| we then present results 
of detailed numerical simulations of BBN in presence of 
matter-antimatter domains. New and stringent limits 



on a ntim atter domains are derived in Sec. VI, whereas 
Sec. VII is devoted to discussion and conclusions. The 



appendices discuss the structure of the actual annihila- 
tion region (App. |A|) and some aspects of the numerical 
treatment of the problem (App. 



II. PRELIMINARIES 

In this work, we consider cosmological models which 
contain different amounts of antimatter distributed in do- 
mains of various sizes r\. In such scenarios, the Universe 
may be envisioned as a distribution of matter with em- 
bedded domains of antimatter, as is schematically shown 
in Fig. [j]. Initially, the matter densities nb(r\r > r&) in 
the matter region and antimatter densities ni{r\r < ta) 
in the antimatter domains with radius ta are assumed 
to be equal throughout the Universe. The average net 
baryon density is thus given by 

n nct = fj aet n-y = n b (r\r > r A )(l - h) - n h (r\r < r A )f h , 

(1) 
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FIG. 1. Sketch of the assumed initial conditions: Antimat- 
ter domains are embedded in a background of matter. 



where n 7 is the average photon number density, fj ne t the 
average net baryon-to-photon ratio, and the filling factor 
fl is defined as the fraction of the volume of the Uni- 
verse filled with antimatter. Here and in the following 
a bar over some quantity denotes the horizon average 
of that quantity. Since baryo-symmetric models will not 
leave baryon number after the completion of annihilation, 
models with an excess of baryon number will be consid- 
ered, i.e. fi < 0.5. Further, we define the antimatter-to- 
matter ratio 



R A = N-JN b 



(2) 



as the ratio of antibaryon number Ni = J y nj(r)d 3 r to 
baryon number N b = J y rib(r)d 3 r. 

It is convenient to express the length scales in our prob- 
lem in comoving units. The length of, e.g., an antimatter 
region at some cosmic time t, or equivalently temperature 
T, may be related to the length it had at a fixed temper- 
ature T , which we choose to be 100 GeV. The physical 
size l(T) of that region in terms of the comoving size 
Zioo = /(100 GeV) is thus given by 



l(T) = I 



100 



V ^ioo / 
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where R(T) is the cosmic scale factor at an epoch with 
temperature T, and we define i?ioo = i?(100GeV) 
I . The time evolution of the scale factor may be de- 
rived from the conservation of entropy, S oc g* s T 3 R 3 = 
const. Thus the scale factor evolves as R oc <?* s T _1 , 
where g* s is the number of relativistic degrees of free- 
dom contributing to the entropy density of the Universe. 
At the electroweak phase transition (T w 100 GeV), 
g* s ~ 100. For definiteness, we will assume that 
g* s (100 GeV) = 100, which allows us to calculate the 
physical length I as a function of the comoving length 



at 100 GeV, Zioo- The comoving scale Zioo thus corre- 
sponds, for example, to a physical length at a temper- 
ature of 1 MeV of f(lMeV) = Zioo x lO^lOO/lOJS) 1 / 3 , 
and to 1{T) = ;i o(100GeV/T)(100/3.909) 1 / 3 , for any 
epoch subsequent to e^-annihilation. 

The time evolution of the densities rii(r,t) of the nu- 
cleon or light nuclei species i (and their antiparticles) 
is governed by three mechanisms, namely diffusion and 
hydrodynamic processes, annihilation, and nuclear reac- 
tions, 



~dt 



dt 



diff /hydro 



dt 



dt 



(4) 



We discuss diffusive and hydrodynamic processes in 
Sec. Ill, whereas baryon- antibaryon annihilations will be 
discussed in Sec. IV. The nuclear reaction network has 
been widely discussed in the literature (e.g. |l^,|l^]) and 
will thus not be described here. Some aspects of the nu- 
merical treatment may be found in App. |b| 

We will frequently use variables Aj(r), which are de- 
fined in terms of the number density of species i at co- 
ordinate r and the average net baryon number density 



Ai(r) 



m(r) 



(5) 



We also will use the quantity S(r) which denotes the rel- 
ative cosmic temperature variation at r, i.e. 



5(r) = 



T{r) - T 
T 



(6) 



where T is an appropriately defined cosmic average tem- 
perature and T(r) is the local temperature. 



III. DIFFUSIVE AND HYDRODYNAMIC 
PROCESSES 

A. Pressure Equilibrium 

Our initial conditions are such that matter and anti- 
matter regions exist in pressure equilibrium with each 
other at uniform cosmic temperature. As long as the 
transport of baryon number over the boundaries from 
one region into the other is not efficient, matter and 
antimatter are kept in separate regions. The photon 
and lepton densities are homogeneous and the temper- 
ature is the same throughout the Universe. Inhomo- 
geneities in the total baryonic density, which is defined 
as the baryon- plus antibaryon- density at position r, 
l n & ot ( r )l = n b( r ) + n l( r )i ma y arise only when transport 
of baryon number over the domain boundaries occurs and 
annihilation proceeds. Subsequently, the baryon and an- 
tibaryon densities close to the boundary decrease, lead- 
ing to a decrease in the (anti-) baryon pressure in the 
annihilation region. This baryonic underpressure is then 
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FIG. 2. Snapshot of an isobaric density fluctuation. The 
full line represents the total baryonic (i.e. baryonic- plus an- 
tibaryonic-) overdensity A(r), the dotted line the deviation 
from the average cosmic temperature S(r). The boundary 
between the matter- and the antimatter- domain is approxi- 
mately at rioo ~ 10 -6 , where rioo is the radius of the spherical 
simulation volume. 



compensated for by a slight adiabatic compression of the 
region and thus an increase of the radiation pressure, 
as to reestablish pressure equilibrium between the re- 
gion and its environment. Fluctuations which have come 
into pressure equilibrium will be termed 'isobaric' fluc- 
tuations [|o). Note that the time scale for reestablishing 
pressure equilibrium is by far the shortest of all time 
scales in our problem, such that the assumption of at- 
taining pressure equilibrium instantaneously is justified. 
An isobaric fluctuation after some annihilation at the do- 
main boundaries has occurred is shown in Fig. [|. 

At late times, and for fluctuations where the photon 
mean free path after e annihilation becomes large com- 
pared to the spatial scale of the fluctuations (typically at 
T s» 20 keV), temperature gradients between the fluctua- 
tions cannot be maintained any more and the assumption 
of pressure equilibrium breaks down. In this regime the 
density inhomogeneities are dissipated by hydrodynamic 
expansion of the cosmic fluid (see below). 



B. Baryon Diffusion 

Transport of (anti-) baryon species i over the domain 
boundaries may only be accomplished by diffusion, which 
is described by 



dru 

at 



(7) 



where D% is the relevant diffusion constant. Equa- 
tion (fjj) may be written in comoving radial coordinates 
(see Ref. pi) 



dAj(T, rioo) 

at 



dr 
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-Ai(T,r W0 ) 
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where we have used the notation introduced in the pre- 
vious section. 

The diffusion constant Dik for baryon species i due to 
scattering off some species k with cross section dik and 
number density nk is approximately given by the product 
of thermal baryon velocity w and baryon mean free path 
lik of the particle under consideration, 



D 



ik 



1 / 11 

^Vil lk = -Vi 

3 6 (JikUk 



(9) 



Some relevant diffusion constants and their cosmologi- 
cal importance may be found in Ref. ]2C| ], The effective 
baryon diffusion constant of nucleus i in the plasma due 
to scattering off different species k is given by 
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E 



i 

~Dik 



(10) 



The diffusion length of a species is defined as the rms 
distance traveled during time t. Written in comoving 
coordinates, one finds (see Ref. pl|j2C| ]) 
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Three different temperature regimes with respect to the 
diffusion of baryons may be distinguished. At early 
times, prior to the annihilation of thermal electron- 
positron pairs, the proton diffusion length is short due to 
their electromagnetic interaction with the ambient pairs. 
Neutron diffusion is controlled by the much weaker mag- 
netic moment scattering on electron-positron pairs and 
thus the diffusion length is longer. As long as the tem- 
perature is higher than w 1 MeV, neutrons and protons 
are however constantly interconverted by the fast weak 
interactions. During this epoch, baryons may thus diffuse 
during the time they spend as neutrons pjj. After the 
weak interactions freeze out at T ss 1 MeV, neutron-to- 
proton interconversion ends, but the neutrons continue 
to diffuse. At T<i He « 80 keV, virtually all free neutrons 
are bound into 4 He nuclei. All baryons and antibaryons 
exist now in the form of charged nuclei and antinuclei. 

Proton and charged nuclei diffusion is limited by 
Coulomb scattering off electrons and positrons from the 
time of weak freeze out down to a temperature T ~ 
40 keV. The Coulomb cross section for the light nuclei 
is proportional to the square of the nuclear charge Zf 
and the thermal velocity to y/l/Ai, where Ai is the mass 
number of the nucleus under consideration. This leads 
to a suppression factor of the diffusivity of nuclei relative 
to that for protons of (l/Zfy/%) @. 

When the pair density decreases at temperatures lower 
than T 40 keV, protons cease to diffuse as individ- 
ual particles. Rather, a proton-electron system diffuses 
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FIG. 3. Diffusion lengths dioo of (anti-) neutrons (dashed 
lines) and (anti-) protons (full lines) as a function of cos- 
mic temperature, measured on the comoving scale fixed at 
T = f 00 GeV. Here a baryon-to-photon ratio of g = 4 x 10~ 10 
has been assumed. 

together in order to maintain electric charge neutrality 
and consequently the larger electron photon cross sec- 
tion dominates the proton diffusion constant J2(J. Diffu- 
sivity of nuclei at these temperatures is expected to be 
suppressed by a factor 1/Zi relative to that of protons. 
The proton and electron diffusion lengths are displayed 
in Fig. |3| for the temperature range of interest. 



C. Heat Diffusion and Hydrodynamic Expansion 

Heat transport between isobaric high and low density 
regions, in particular between regions with high and low 
total baryonic pressure and concomitant low and high 
radiation pressure, may be accomplished by diffusing or 
free streaming neutrinos or by diffusing photons. Note 
that the effect of such heat transport is the decrease of 
(anti-) baryonic density in overdense regions, and the in- 
crease of (anti-) baryonic density in underdense regions. 
We will be interested in the evolution of inhomogeneities 
which are generated by annihilations at relatively low 
temperatures T < a few MeV, such that neutrino heat 
conduction is typically inefficient pOfl . In contrast, heat 
transport via photons may be efficient towards the end 
of the annihilation, depending on the length scale of 
the temperature fluctuations. During the time of an- 
nihilation, the comoving photon mean free path 
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l 100 



R- 



<JTn e ± 



(12) 



increases enormously since it is inversely proportional 
to the total number density of electrons and positrons, 
n e ± = n e + + n e - . At temperatures below T < 30 keV 
essentially all pairs have annihilated and n e ± is domi- 
nated by net electron number densities required to main- 
tain charge neutrality. The increased photon mean free 



path may then affect the dissipation of fluctuations in 
the baryon density [f23| . As long as the photon mean 
free path is still shorter than the scale of the fluctuation, 
heat transport is described by the diffusion equation for 
photons, which is identical to Eq. (fjj), but with A^(r) 
replaced by the temperature fluctuations S(r). The dif- 
fusion constant is now given approximately by 
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(13) 



with g t the statistical weight of the heat transporting par- 
ticles (gt = 2 for photons) and g* the statistical weight 
of the relativistic particles still coupled to the plasma 
(<?* = gt after e ± annihilation, since neutrinos are decou- 
pled). 

When the photon mean free path becomes larger than 
the scale of fluctuations, free-streaming photons will keep 
high- and low- density regions isothermal with baryonic 
pressure gradients remaining. In this regime, dissipa- 
tion of inhomogeneities proceeds via expansion of high 
density regions towards low density regions and the con- 
comitant transport of material towards the annihilation 
region. The motion of the charged particles, protons 
and light elements, is impeded by the Thomson drag 
force, which acts on the electrons dragged along by the 
charged nuclei [|4j. Balance between pressure forces 
and the Thomson drag force yields a terminal velocity 
v = dr wa /dt @, 



1 dP 



(14) 



where dP/drioo is radial (anti-) baryonic pressure gradi- 
ent, rioo radial coordinate as measured on the comoving 
scale, e 7 photon energy density, and n e = n e - — n e + net 
electron density. One finds for the pressure exerted by 
baryons and electrons below T ss 30 keV 



P » Tn^ (J2 A, + (n; air 2 + (]T Z { A,) ' 



(15) 



with the sum running over all nuclei i with nuclear charge 
Zi . Note that expression ( Jl5| ) quickly reduces to the pres- 
sure exerted by an ideal gas when the reduced pair 
/n nct becomes negligible compared 



density n* a 
to J2i ZiAi 
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IV. MATTER-ANTIMATTER ANNIHILATION 

A. Annihilation Reactions and Cross Sections 

The dominant process in nucleon-antinucleon interac- 
tion is direct annihilation into pions, 
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FIG. 4. Interaction rates for pions with leptons (dotted 
lines) and nucleons (full lines) as a function of cosmic tem- 
perature. Here a baryon-to-photon ratio r\ = 4 x 10~ 10 has 
been assumed. For comparison, the pion decay rates are also 
shown (dashed lines). 

Electromagnetic annihilation (p+p — > 7+7) is suppressed 
by a factor of (m e /m p ) 2 « 3 x 10~ 7 . Annihilation via 
the bound state of protonium is also possible, but the 
cross section is smaller by (m e /m p ) 3 ^ 2 w 10 -5 compared 
to direct annihilation. 

The charged pions either decay with a lifetime of r n ± = 
2.6 x 10~ 8 s directly into leptons 



7T — ► fT + v p 

L 



e + i/ e -+- i/, 
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e + z/ e + za. 



or may be transformed into tt° via charge-exchange 
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or weak interactions 
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The neutral pions subsequently decay into photons with 
T n a = 8.4 x 1CP 17 s. The rates for the three channels, de- 
cay, charge exchange, and weak interactions are shown in 
Fig. [|. Below a temperature of a few MeV, decay dom- 
inates the loss of charged pions unless the local baryon- 
to-photon ratio is well in excess of 77 ~ 1CP 10 . Never- 
theless, even at low rj some pions may charge-exchange 
on nucleons; possible consequences thereof will be dis- 
cussed below. Weak interactions with the ambient lep- 
tons do not significantly contribute to the pion interac- 
tion rates in the temperature range relevant for our work 



TABLE I. Probabilities for the production of secondary 
nuclei and nucleons in p 4 He annihilations, derived from the 
branching ratios given in Ref. J25[] . 
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0.51 


0.28 


0.13 


0.43 


0.21 



(T < 20 MeV). Neutral pions never have a chance to in- 
teract with either leptons or nucleons, due to their rapid 
decay. 

Annihilation of antinucleons on light nuclei, TV, pro- 
duces a wealth of secondary particles, 



p/n + N — > N ,p,n,ir 



(21) 



In Table | we give probabilities for the production of 
secondary nuclei N in the 4 He annihilation process fol- 
lowing Ref. |p5f . Production of secondary nuclei during 
antinucleon annihilations on the other light nuclei, and 
anti- 4 He annihilation on 4 He, are relatively unimportant 
for this study. Due to the huge difference in the abun- 
dances of 4 He on the one hand and the other nuclei on 
the other hand, the destruction of only a minute frac- 
tion of 4 He may already have significant impact on the 
abundances of the other primordial elements. Disrup- 
tion of the other elements will occur with much smaller 
probabilities due to the smaller abundances of these iso- 
topes, and the secondaries of these processes will never 
contribute significantly to the respective abundances. In 
our numerical calculations wc assumed that disruption of 
all elements but 4 He results in free nuclei only. 

We are interested in annihilations of antimatter do- 
mains between shortly before the epoch of weak freeze 
out and recombination such that the annihilation cross 
sections for thermal nucleons with kinetic energies be- 
tween a few MeV and about 10 -7 MeV are needed. Ex- 
perimental data are available only down to an incident 
momentum of about 30-40 MeV, corresponding to ki- 
netic energies of about 1 MeV. Therefore, we have to 
extrapolate the experimental data with the help of ex- 
isting theoretical calculations for the cross sections down 
to the relevant energy range. At such low energies, the 
Coulomb forces between charged particles become impor- 
tant, thus systems with Coulomb interactions like pp and 
pN, and without, like fin, rip, and fiN, have to be treated 
separately. 

The product of annihilation cross section cr ann and rel- 
ative velocity v in systems with at least one neutral par- 
ticle is known to be approximately constant at low ener- 
gies [ p6[ . Experimental values of cr ann w = 40±3mbc and 
32 ± 5 mb c were obtained at center of mass momenta of 
22 MeV/c and 43 MeV/c, respectively Q. In our cal- 
culations, we used a constant value of a ann v = 40mbc. 

In systems with Coulomb interactions, such as the pp 
or the pN system, the behavior of the low-energy an- 
nihilation cross section is drastically modified due to 
Coulomb attraction. Indeed, the charged particle low 



G 



energy annihilation cross section is found to be inversely 
proportional to the square of the incident momentum 
and therefore the reaction rate is formally divergent at 
zero energy. (This divergence is of course removed when 
matter and antimatter reach atomic form at low temper- 
atures.) Again, there are no experimental data below 
about 1 MeV kinetic energy. The available experimen- 
tal data at higher energies is however well reproduced 
by phenomenological calculations pOpl| ]. The results of 
these calculations depend on the mass and the charge 
of the nucleon/nucleus under consideration as well as on 
some phenomenological parameters, which may be deter- 
mined experimentally. 

For energies below about 10 _1 MeV, which are mostly 
relevant for our study, the thermally averaged pp annihi- 
lation rate in a plasma of temperature T may be approx- 
imated by 



< cr ann v > w f Pp mb c 



MeV 



T 



(22) 



where fp p = 32 is a numerical constant. The equiva- 
lently defined numerical constants for the pN annihila- 
tion channels are of similar magnitude, /^h ~ /p 3 H ~ 16 
and /p3 He ~ /p 4 He ~ 20. These were obtained using the 
results of the phenomenological calculations as given in 



Refs. [26, 
Refs. [28 1 



based on the experimental data as given in 



B. Impact of Secondaries in Nucleon-Antinucleon 
Annihilations on BBN 

It is of interest if annihilation-generated photons and 
pions, or their decay products, may alter the abundance 
yields, either through their effect on weak freeze out or 
by, for example, photodisintegration or charge exchange 
reactions. In a single annihilation event, about 5-6 pions 
with momenta ranging from tens to hundreds of MeV are 
produced. Reno & Seckel j3^[ showed that prior to 
annihilation the thermalisation time scale for charged pi- 
ons is always shorter than the hadronic interaction time. 
At later times, the charged pions have no chance to in- 
teract due to their short lifetime. We may thus use the 
thermally averaged charge exchange cross sections (cf. 
Eq. |l^) given in Ref. plj to calculate the ratio between 
a typical charge exchange interaction time and t^± 



7"cex 
7V ± 



0.01 



?7local 

4 x 10- 



— V 

MeV / 



(23) 



Charged pions may thus only charge exchange, if ?7iocai 
10~ 10 , and the temperature is not much lower than 
1 MeV. Due to their electromagnetic interaction, charged 
pions remain mainly confined to the annihilation region 
as long as pairs are still abundant. Within that region, 
?7iocal is typically much lower than the (anti-) baryon- 
to-photon ratios anywhere else due to prior matter- 



antimatter annihilation (cf. App. ^). At lower tempera- 
tures, when the pions may easily move within the primor- 
dial fluid, charge exchange reactions are negligible due to 
the small nuclear densities, as is apparent from Eq. 
We therefore observe only negligible impact of charge ex- 
change reactions on the BBN abundance yields in our nu- 
merical simulations, even for early matter-antimatter an- 
nihilation and large Ra (implying large 7yi OC ai in matter- 
and antimatter- domains). 

The leptonic secondaries ^ il ,e ± and v do not mod- 
ify the details of weak freeze out, unless the number of 
annihilations per photon is extremely large. As long as 
this number is not approaching unity, annihilation gen- 
erated v e 's have negligible effect on the n/p-ratio, since 
their number density is orders of magnitude smaller than 
that of the thermal f e 's, which govern the weak equilib- 
rium. The same holds for electrons and positrons pro- 
duced in /x-decay which are quickly thcrmalized by elec- 
tromagnetic interactions. 

In each annihilation event about half of the total en- 
ergy is released in form of electromagnetic energy. Impor- 
tant impact on the BBN abundances may result through 
these annihilation-generated 7-rays and e cascading on 
the background photons (and on pairs before e anni- 
hilation) via pair production and inverse Compton scat- 
tering on a time scale rapid compared to the time scale 
for photodisintegration of nuclei [p3|-p5[. After cosmic 



e annihilation, the cascade only terminates when in- 
dividual photons do not have enough energy to further 
pair-produce on background photons. For temperatures 
T > 5 keV, the energy of 7-rays below the threshold 
for e ± -production does not suffice for the photodisinte- 
gration of nuclei. If annihilations occur below 5 keV, the 
light nuclei gradually become subject to photodisintegra- 
tion, according to their binding energy. 

Destruction of 2 H, 3 H and 3 He by photodisintegration 
is thus possible for temperatures below T < TL 2 h « 



5 keV and T < TL,3 He w T" 7 3 H w 2 keV, respectively. 
This destruction is nevertheless subdominant to the pro- 
duction of these isotopes by photodisintegration of 4 He, 
possible at lower temperatures (T 7 4 Hc < 0.4 keV), sim- 
ply due to the larger abundance of 4 He. Thus photo- 
disintegration of 2 H, 3 H and 3 He may only be impor- 
tant if annihilations take place in the temperature range 
0.4 keV < T < 5 keV. However, direct production of 2 H, 
3 H and 3 He via annihilations on 4 He dominates destruc- 



tion by photodisintegration (see Sec. IV C below). The 
photodestruction factor for, e.g. 2 H, may be roughly es- 
timated from the photodestruction factor for 4 He given 
in Ref. |3j|. Since the cross section for the competing 
process, i.e. Bethe-Heitler pair production on protons, 
does only slowly vary with temperature, these destruc- 
tion factors may be used in the relevant temperature 
range, 5keV < T ann < 0.4 keV, but have to be scaled 
by the target abundances. Thus one estimates a fraction 
of about 0.1 ( 4 He/ 2 H) » 10~ 4 2 H nuclei per GeV electro- 
magnetic energy injected to be photodisintegrated. Di- 
rect annihilations will create about 0.2( 4 He/p) « 10~ 2 
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2 H nuclei per annihilation. Even though a fraction of 
those will thermalize within the annihilation region, and 
thus be subsequently annihilated, it is well justified to 
neglect photodisintegration of lighter isotopes. Photodis- 
integration of nuclei lighter than 4 He was thus not taken 
into account in our simulations. 

For the 4 He photodisintegration yields we used the re- 
sults given in Protheroe et al. Pqj . Production of 3 He 
and 3 H exceeds the 2 H yield by a factor of ten, thus 
typically producing a large 3 He/ 2 H ratio. The photodis- 
integration yield for 3 He is peaked at a temperature of 
T w 70 eV and becomes significantly smaller at lower 
temperatures. When photodisintegration of 4 He occurs 
it creates initially energetic 3 He and 3 H nuclei. The im- 
portance of the energetic 3 H and 3 He nuclei resulting 
from photodisintegration is twofold; they do not only di- 
rectly increase the 3 He abundance, but may also lead to 
production of 6 Li via 3 H/ 3 H + 4 Hc — > 6 Li + n/p, as 
was recently stressed by one of us [[36). We take th e 6 Li 
yields as given in Ref. |36| (for more details, cf. Sec. $VC\ 
below). 

Note that all photodisintegration yields have been 
calculated using the generic 7-ray spectrum given in 
Protheroe et al. |35|. This procedure is only adequate 
as long as the energy of the injected photons, E-y « 
200 MeV, is beyond the threshold for e ± -pair produc- 
tion of the injected 7's on the background photons, 
E c (z) ~ 4.7 x 10 7 (1 + z)- 1 MeV, i.e. for z > 2 x 10 5 , 
corresponding to T ann > 50 eV. The results for later anni- 
hilation, i.e for antimatter domains on scales larger than 
r A°° It 10 3 cm thus have to be interpreted with some 
care. In order to be able to give conservative limits, we 
have done simulations where photodisintegration was ig- 
nored which gave weaker limits by a factor of a few (see 
Fig. lift. 



average, the secondary nuclei gain a kinetic energy Eq 
of a few tens of MeV [^5| . Their transport is then ini- 
tially described by free-streaming until their energy has 
decreased to thermal energies through interactions with 
the plasma, after which transport has to proceed via ther- 
mal diffusion. The dominant energy loss mechanisms for 
energetic charged nuclei in a plasma with kinetic energy 
below 1 GeV are plasmon excitations and Coulomb scat- 
terings. Note that these processes have negligible impact 
on the direction of the momentum of the energetic nuclei 
such that free-streaming is a good approximation [ p7[ . 
The distance covered until the kinetic energy of the par- 
ticles has decreased to the thermal energy of the plasma 
defines the stopping length, 



^stop 



jstop 



dx 



dx 
dE' 



-dE 



(24) 



Provided ^ stop is larger than the size of the annihilation 
region, all nuclei which become thermalized in a mat- 
ter domain (typically about 1/2) have a good chance to 
survive. If we calculate the energy loss per distance fol- 
lowing Ref. m we find for the stopping length measured 
in our comoving coordinates 



jstop 

noo 



55 cm 



/keV 
4 x 10" 10 



50 MeV 



(25) 



where Z is the charge of the energetic nucleus. In an 
analogous manner we may calculate a stopping time 



'stop 



-^thermal 



1 da; 
v{xjdE 



dE 



(26) 



C. Impact of Secondaries in Antinucleon-Nucleus 
Annihilations on BBN 

While we expect secondaries of nucleon-antinucleon 
annihilations to have only significant effect at lower tem- 
peratures, energetic nuclei arising in antinucleon-nucleus 
annihilations may substantially modify the light element 
abundances for temperatures as high as T < T4 Ho w 
80 keV. This may occur through direct production of light 
isotopes in antinucleon annihilations on 4 He as well as 
through possible subsequent non-thermal fusion of these 
energetic light isotopes on 4 He. Since the 4 He abun- 
dance exceeds the abundances of the other isotopes by 
orders of magnitude, p 4 He is the dominant antinucleon- 
nucleus annihilation process. The relative probabilities 
for the production of the various secondary nuclei arising 
in 4 He disruption are given in Table [|. The secondary 
nuclei are produced within the annihilation region and 
may thus themselves be subject to annihilation, unless 
they are able to escape from the annihilation region. On 



needed to slow down a particle to thermal energies. Eval- 
uation of the integral yields for charged particles 



stop 



3.25 s 



/keVV 
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3/2 



E, 



3/2 

thermal 



\ T J 
4 x 10 

?7local 



(50 MeV) 3 / 2 



10 



Z~ 



(27) 



Neutrons lose their energy through nuclear scatterings. 
In contrast to the charged particle interactions discussed 
above, the deflection angle in a nuclear scattering event 
may be large, such that the use of Eq. ( p4] ) is inappro- 
priate since it relies on the free-streaming assumption. 
Rather, the distance covered by the neutrons is described 
by a random walk. The stopping time is nevertheless de- 
scribed by Eq. (|2^), since the energy loss does not depend 
on the direction of the motion. The energy loss per unit 
distance for neutrons may be estimated via 



dE 

— = {- h\j)E(j np n p , 



(28) 



8 



where / is an approximate average fractional energy loss 
in each elastic neutron-proton scattering event. If we as- 
sume a simple power law for the neutron-proton cross sec- 
tion a np w 10 3 mb(S/10MeV) 115 (cf. Fig. 1 in Ref. |§) 
and an energy loss of 80 % in each scattering event, we 
find 



stop 



1.56 x 1(T s 



keV 



4 x 10- 

Tfrocal 



p0.65 
^0 



p0.65 
^thr 



(10 MeV) - 65 



(29) 



We may now compare r st0 p for neutrons and pro- 
tons with a typical 4 He spallation time scale, T sp = 
((<7 sp i>)n4 Hc ) . We find that only about a fraction of 
10~ 3 of all energetic protons may spallate additional 4 He. 
Since direct production of energetic protons and light nu- 
clei in p 4 He disruption is of similar magnitude, this effect 
may safely be neglected. For the energetic neutrons, we 
find that about 30 % may spallate 4 He. Since we ob- 
tain energetic neutrons in about half of the annihilation 
events, additional 3 H or 3 He will be produced in about 
a tenth of the p 4 He annihilations. This is not signifi- 
cant compared to the secondary 3 H and 3 He produced 
directly in the p He annihilations with a probability of 
about 60% and will therefore be neglected in the numer- 
ical treatment. 

The annihilation generated energetic light nuclei may 
however be important as a source for very rare light el- 
ements such as 6 Li via the fusion reactions 3 H + 4 He 
-> 6 Li + n and 3 He + 4 He -> 6 Li + p §|. Using a 
value of 35 mb for the fusion cross section, the thresh- 
old energies of E th = 4.80 MeV and E th = 4.03 MeV, 
respectively, and the energy distribution for the nonther- 
mal mass three nuclei as given in |25J, we find 



3 H 4 Hc-*n 6 Li/ 



2 x 10" 



and 



(P, 



He 4 He^p 6 Li 



) w 5 x 10" 



(30) 



(31) 



for the probabilities to produce 6 Li from energetic mass 
three nuclei. The calculation is done similarly to the ones 
in Refs. p8| , |3"6[ ] , where energy loss of energetic nuclei ac- 
cording to Eq. ( pjj ) has been taken into account. The 
number of 6 Li nuclei produced per antiproton annihila- 
tion is thus 



X 



°"p 4 Ho 

1.8 x 10~ 8 



™ 4 Hc 



0.25 J ' 



where and P3n e are the probabilities to create 3 H or 
3 He in a p 4 He annihilation event (see Tab. m. A simple 
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FIG. 5. Spatial dimension of the antimatter regions and 
corresponding annihilation temperature for an antimatter 
fraction of -Ra = 0.1. Also given is the baryonic mass in units 
of the solar mass M© contained in the antimatter region and 
the annihilation redshift. Note that annihilation is stalled in 
the temperature range T4 Hc > T > 5keV, as is explained in 
the text. 

estimate for the total synthesized 6 Li/H abundance (ex- 
cluding production via 4 He photodisintegration) is thus 



HA 
H 



1.8 x 10~ 9 



nb J 



Ra 
0.1 



0.25 



(33) 



for i?A *C 1 and where it is understood that only that 
fraction of antimatter has to be inserted in Eq. ( p3[ ) which 
has not annihilated by temperature T4 He ps 80 keV. This 
is many orders of magnitude higher than the standard 
BBN value, n6 Li /n p = O(10~ 13 ) and will therefore pro- 
vide very stringent limits in some areas of the parameter 
space, as will be discussed in Sec. [v|. 



V. BBN WITH MATTER-ANTIMATTER 
DOMAINS 

After having discussed the different dissipation mech- 
anisms of antimatter domains in the early Universe as 
well as the annihilation reactions and the possible im- 
pact of annihilation generated secondaries on BBN, we 
are now in a position to put all this together in order 
to examine the influence of annihilating antimatter do- 
mains in the early universe on the BBN light element 
abundance yields. Clearly, such scenarios involve such a 
multitude of nuclear reactions and hydrodynamic dissipa- 
tion processes that obtaining fairly accurate predictions 
for the BBN yields requires numerical simulation. We 
have therefore substantially modified the inhomogeneous 
BBN code by Jedamzik, Fuller & Mathews pi] , origi- 
nally including nuclear reactions, baryon diffusion, and 
fluctuation dissipation by photon and neutrino induced 
processes, as to also include nuclear reactions between an- 
timatter, matter-antimatter annihilation reactions, free- 
streaming of secondary nuclei produced in annihilations, 
the non-thermal fusion reactions of secondaries, as well 
as photodisintegration of 4 He through annihilation gen- 
erated cascade 7-rays. Some processes are not included 
in our simulations, due to their marginal impact on BBN 
as outlined in the last section. For more details on the 
procedure of the numerical simulation the reader is also 
referred to App. H. 
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FIG. 6. Results of BBN calculations in a Universe with antimatter domains. Shown are the 4 He mass fraction, the abundances 
of 2 H, 3 He, 6 Li and 7 Li relative to hydrogen and the ratio of 3 He over 2 H as a function of antimatter domain radius, r^°, for 
several low values of the antimatter-to-matter ratio, Ra (see legend). 



A detailed analytic and numerical analysis of the ac- 
tual structure of the annihilation region, i.e. the region 
at the domain boundaries where the bulk of annihilations 
occur, is presented in App. |A|. Our simulations do not 
accomplish to resolve the physical width, Z ann , of this re- 
gion due to the increasingly large ratio, Ra/ l&nm and the 
required extreme dynamic range to resolve both scales. 
However, we believe, as we argue in App. |A|, that this 
"flaw" of our simulations has only little impact on the 
accuracy of our results. This is also confirmed by the 
relative independence of our results on the total number 
of zones employed in the simulations (see App. |Bj). 

The relationship between antimatter domain size, r\ 00 , 
and approximate annihilation temperature T and red- 
shift z of a domain is shown in Fig. ^| and is determined 
by neutron diffusion at early times and hydrodynamic 
expansion at late times. Since dissipation by neutron 
diffusion for temperatures T > T4 Hc sa 80keV is rela- 
tively more efficient than dissipation by hydrodynamic 
expansion at somewhat lower temperatures, antimatter 
domain annihilation does typically not occur in the tem- 



perature regime between T4 He and 5 keV. This implies 
that "injection" of antimatter and annihilation between 
the middle and the end of the BBN freeze-out process, as 
envisioned in the scenarios of Ref . |l6[ , may not operate 
in scenarios with a matter-antimatter domain structure 
in the early universe. 

We present the detailed numerical results of our study 
in Figs. |^ and ^. Shown are the abundances of the re- 
spective elements and the 3 He/ 2 H ratio for a number of 
matter-antimatter ratios Ra as a function of the typ- 
ical size of the antimatter regions r A 00 . In the subse- 
quent discussion of our results we will distinguish be- 
tween three different limiting cases, according to the 
segregation scale of matter and antimatter, or equiva- 
lently, the approximate matter-antimatter annihilation 
time. The fractional contribution of baryons to the crit- 
ical density of the Universe, Qb, was kept at a constant 
value of SI5 = 0.0125/ijqq in all simulations, where h wo 
parameterizes the value of the Hubble parameter Hq to- 
day, ftioo = -Ho/(100 km s -1 Mpc -1 ). The corresponding- 
value of the SBBN parameter is rj = 3.4 x 10~ 10 . 
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A. Annihilation Before Weak Freeze Out 

The key parameter which determines the primordially 
synthesised amount of 4 He is the n/p-ratio at a temper- 
ature of T4 Hc ss 80 keV. Early annihilation of antimatter 
proceeds mainly via neutrons diffusing towards the an- 
timatter domains and being annihilated at the domain 
boundaries whereas antineutrons diffusing towards the 
matter domains may annihilate on both, neutrons and 
protons. The net effect of both processes is the prefer- 
ential annihilation of antinucleons on neutrons, perturb- 
ing the n/p-ratio towards smaller values as detailed in 
Rcf. jn|. (One finds this effect often to be even more pro- 
nounced since after some initial annihilation of protons 
close to the domain boundary the annihilation region is 
completely void of protons.) If these perturbations in the 
n/p-ratio persist down to T4 He , a significant reduction in 
the synthesized 4 He mass fraction may result. A poten- 
tial effect in the reverse direction, increase of the n/p- 
ratio by preferential pion-nucleon charge exchange reac- 
tions on protons (cf. Eq. |l9| ) , is subdominant as discussed 
in Sec. IV B. It is of interest at which temperatures these 



possibly large perturbations in the n/p-ratio may still 



be reset by proton-neutron interconversion via weak in- 
teractions governed by the rate r wca k « G| ermi T 5 . In 
the upper two panels of Fig. || we show the n/p-ratio as 
a function of temperature for comparatively early anti- 
matter domain annihilation. In panel (a), the antimatter 
fraction was chosen to be Ra = 0.5 and the length scale 
of the antimatter regions to be r A 00 = 0.018 cm, corre- 
sponding to an approximate annihilation temperature of 
about 5 MeV. The n/p-ratio for this parameter combi- 
nation is observed to be virtually indistinguishable from 
the n/p-ratio in a standard BBN (SBBN) scenario. Thus 
the final 4 He mass fraction (dashed line) coincides with 
the SBBN value of Y p w 0.24. Only when the matter-to- 
antimatter ratio rises to values of order unity (panel b), 
the n/p-ratio in the presence of antimatter annihilations 
deviates significantly from the corresponding quantity in 
a SBBN scenario (shown by the dotted line). However, 
after the annihilation has completed, the weak interac- 
tions are still rapid enough to reestablish weak equilib- 
rium and thus the final 4 He mass fraction and the other 
light element abundances emerge unaffected. In Figs. ^ 
and ^ one observes that significant impact on the synthe- 
sized 4 He mass fraction occurs only for r A 00 larger than 
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5 x 10 2 cm to 10 'cm, depending on Ra- 



B. Annihilation After Weak Freeze Out, But Before 
4 He Synthesis 



The situation changes when annihilation occurs dur- 
ing or after weak freeze out, when neutron-proton inter- 
conversion ceases to be efficient. Annihilation continues 
to proceed mainly via neutrons and antineutrons, since 
proton diffusion is still hindered by the abundant e ± 
pairs. However, neutrons and antineutrons which have 
diffused out of their respective regions may now not be 
reproduced anymore. Antimatter regions of typical size 
larger than the neutron diffusion length at the cosmo- 
logical epoch of weak freeze out thus provide a very ef- 
ficient sink for neutrons. The n/p-ratio is strongly af- 
fected in such models as is apparent from the lower two 
panels of Fig. |[ In a model with, e.g., Ra = 0.1 and 
r A°° = 0-55 cm J annihilation proceeds at a temperature 
of m 0.8 MeV (panel c). At this temperature, the weak 
interactions are not rapid enough to reestablish the equi- 
librium n/p-ratio. Thus the final 4 He mass fraction Y v 
will be decreased compared to its SBBN value. This is 
also evident from Figs. ^ and |?] for antimatter domains 
with length scales between ~ 5 x 10 -2 cm and about 6 cm. 
For small antimatter fractions, Ra ^$ 0.1, the other light 
element abundances are comparatively less affected than 



4 He, only for larger values of Ra production of 2 H, 3 He, 
and 7 Li is also strongly suppressed. 

Since the primordial 4 He mass fraction mostly depends 
on the n/p-ratio at T4 Hc , it may be estimated analytically 
whenever n/p|T 4He is known, 
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If one assumes that annihilation occurs instantaneously, 
the n/p-ratio in a scenario with annihilating antimatter 
domains may be estimated by 



{n /n b ) exp [—Ati/T n ] - xR A 
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x exp [-At 2 /r n ] , 



x)Rji 



(34) 



where x is the fraction of antibaryons annihilating on 
neutrons, no and po are the (pre-annihilation) neutron 
and proton densities at T « 0.2 MeV, and nt, is the ac- 
tual baryonic density in the matter region. Neutron de- 
cay is taken into account by the two exponentials, where 
Ati is the time interval between the moment after which 
the neutron fraction is (apart from annihilations) only af- 
fected by neutron decay (T w 0.2 MeV) and the moment 
of annihilation, while At 2 is the time remaining until 
neutrons are incorporated into 4 He at T4 Ho ss 80 keV. 
Thus the two limiting cases between which this estimate 
should hold are identified by At\ « s, At 2 ~ 130 s (an- 
nihilation at T w 0.2 MeV) and At x « 130 s, At 2 ~ s 
(annihilation at T w 80 keV), respectively. Note that 
Eq. ( [34| ) neglects the increase in proton density due to 
neutron decay. The fraction x is well approximated by 
x rs 1, since there are practically no protons present in 
the annihilation region such that most antibaryons an- 
nihilate on neutrons. The results of the above estimate 
agree remarkably well with the numerical results. 

It is apparent from Eq. (HJ) that the n/p-ratio at T4 Hc 
not only depends on the antimatter fraction Ra, but also 
on the time when annihilation of the antimatter domains 
takes place. The reason for this behavior is that the 
number of neutrons annihilated is roughly independent 
of the annihilation time, but for early annihilation this 
number is subtracted from a larger initial number than 
in case of later annihilation. 

The above estimate Eq. (|3) also predicts that it is 
possible to completely avoid^Ie, and in that case also 
2 H, 3 H, 3 He, and 7 Li synthesis, namely if n/p|T 4Hc = 0. 
(For antimatter fractions which yield negative results for 
n/p Eq. ([li]) is obviously not applicable.) Thus there 
is no lower limit to the production of He. An exam- 
ple of such a scenario is shown in panel (d) of Fig. ||. 
The antimatter fraction of Ra = 0.2 exceeds the neutron 
fraction at the time of annihilation (« 0.5 MeV) and thus 
practically all neutrons are annihilated and no light ele- 
ment nucleosynthesis is possible (cf. also the results for 
Ra ^ 0.2 shown in Fig. 0). This is to our knowledge the 
only baryo-asymmetric scenario in which light-clement 
nucleosynthesis is absent [ ju| . 
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C. Annihilation After 4 He Synthesis 

Essentially all free neutrons and antineutrons will be 
bound into 4 He and anti- 4 He nuclei at T4 Hc . At this time, 
the neutron diffusion length is about 6 cm, thus antimat- 
ter domains which are larger than w 6 cm may not an- 
nihilate induced by mixing of matter and antimatter via 
neutron diffusion. Further annihilation is delayed until 
transport of protons and light nuclei over this distance 
is effective. The proton diffusion length does not grow 
to this value until the temperature drops down to a few 
keV. But at this low temperature, the photon mean free 
path has already increased enormously and thus bary- 
onic density gradients in the primordial fluid may not be 
supported any more by opposing temperature gradients. 
Thus the regions far away from the annihilation region, 
which are at high (anti-) baryon density, quickly expand 
towards the baryon depleted and thus low-density annihi- 
lation region and thereby transport matter and antimat- 
ter towards the boundary. The annihilation time is thus 
controlled by the hydrodynamic expansion time scale at 
late times. Only the actual mixing, i.e. the transport 
over the boundary, is still described by baryon diffusion. 

During the course of late time annihilation not only 
nucleon-antinucleon, but also antinucleon-nucleus anni- 
hilations may take place. The elemental abundances pro- 
duced at the BBN epoch may now be substantially mod- 
ified not only by direct annihilations, but also due to the 
effects of the secondary nuclei produced in antinucleon- 
nucleus annihilations. In particular, annihilations on He 
produce 2 H, 3 H, and 3 He nuclei, which, since they are 
energetic, may fuse vi a no n-thermal nuclear reactions to 
form 6 Li (cf. Section IV). Furthermore, the photodis- 
integration of 4 He by energetic photons arising in the 
annihilation process becomes possible. These effects are 
evident from Figs, p] and [?]. Whenever r\ 00 > 6 cm, the 
yields for 2 H and 3 He, as well as for 6 Li show a strong 
increase. The abundance of 7 Li is not much affected by 
late time annihilation, since for Ra not too large there 
is no efficient production channel leading to this isotope, 
and destruction via direct annihilation is insignificant. 
The slightly elevated value for 7 Li/H compared to the 
SBBN result is only due to our initial conditions which 
lead to a higher value of r\ during the BBN epoch for late 
annihilation of high antimatter fractions. ( 7 Li/H is an 
increasing function of r\ for r\ > 3 x 10~ 10 .) 

Since annihilation of antimatter in a scenario where an- 
timatter is distributed in well-defined domains is mainly 
confined to the region close to the matter-antimatter 
boundary, one may speculate that secondary nuclei which 
are produced inside the annihilation region are also an- 
nihilated. But this is not necessarily the case. The sec- 
ondary nuclei gain a kinetic energy of the order of a few 
tens of MeV in the 4 He disruption process pa]. The 
fraction which may survive the annihilation depends on 
the distance the nuclei may free-stream away from the 
boundary before they thermalize. This fraction is ini- 
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FIG. 9. Abundance yields for 2 H and 3 He + 3 H (upper 
row), and 6 Li (lower row) obtained in consideration of differ- 
ent physical effects. Left column: all effects included. Mid- 
dle column: only production by photodisintegration (dashed 
line) and only direct production by annihilation and escape 
from the annihilation region (solid line) taken into account. 
Right column: only direct production by annihilation is taken 
into account, secondaries remain either confined to the an- 
nihilation region (dashed line) or are homogeneously dis- 
tributed throughout the simulation volume (solid line). See 
text for detailed discussion. The antimatter fraction in the 
simulations shown here is -Ra = 10 -3 and the length scale 
r A 00 = 5.5 x 10 2 cm. 



tially small but increases rapidly due to an increase of the 
stopping length for energetic nuclei with cosmic time. (Of 
course, nuclei which travel into the antimatter domains 
and thermalize there will always be annihilated.) In par- 
ticular, according to Eq. ( |25| ) an energetic 3 He nucleus 
of 50 MeV produced via annihilation at T w 3 keV has a 
stopping length of ~ 1.5 cm. This should be compared 



to the typical domain size 



„ioo 



10 cm of domains an- 



nihilating at T ~ 3keV (cf. Fig. [3]), illustrating that the 
synthesized 3 He nuclei will predominately be annihilated 
subsequently. In contrast, when the annihilation occurs 
at T w 60 keV, for domain size r A 00 « 10 3 cm, the stop- 
ping length has increased to«4x 10 3 cm implying that 
most annihilation produced 3 He survives. These trends 
are evident from Figs. ^| and |?] where for the same Ra 
one observes a rapid increase of production of 2 H, 3 Hc, 
and 6 Li with increasing antimatter domain size. The in- 
crease is enhanced also due to the additional production 
of 2 H, 3 He, and 6 Li via photodisintegration of 4 He which 
becomes possible at low temperatures. 

In order to gauge the relative importance of the two 
effects which yield energetic 3 He nuclei, namely direct 
annihilation on 4 He and photodisintegration of 4 He, we 
show in Fig. || results for simulations with Ra = 10~ 3 
and r A 00 = 5.5 x 10 2 cm, in which annihilation occurs 
close to the temperature where the photodisintegration 
yields are maximal (T w 50 eV). Shown are the abun- 
dance yields for 2 H and 3 He (upper row) and for 6 Li 
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(lower row). In panels (a) and (d) all effects, production 
by photodisintegration as well as direct production by 
annihilation and escape from the annihilation region, are 
included. In the middle column, panels (b) and (e), only 
one of the mechanisms is active at a time. The solid line 
shows the results for 2 H, 3 He, and 6 Li when only direct 
production by annihilation is taken into account, while 
the dashed line only considers 4 He photodisintegration 
induced production of these elements. We find that the 
photodisintegration yields for 3 He are larger by nearly a 
factor of 2 than the direct annihilation yields. This is 
not surprising, given the peak photodisintegration yield 
of about 0.1 3 He nuclei per annihilation (cf. f55f) and 
the probability of direct 3 He production in a ^Frle an- 
nihilation weighted by the relative abundance of 4 He to 
protons; P3 He +3 H ( 4 He/p) w 0.05 (cf. Tab. |). For the 
production of 6 Li from energetic 3 He and 3 H nuclei both 
effects arc of the same importance, since the yields for 
6 Li production via energetic nuclei generated by photo- 
disintegration of 4 He and annihilation of 4 He are similar 
over a wide range of redshifts, 2 x 10 4 < z < 4 x 10 5 . 
The remaining two panels, (c) and (f), demonstrate once 
more the importance of the escape of the energetic sec- 
ondaries from the annihilation region. The dashed line 
shows the results of a simulation where all annihilation 
generated nuclei are confined to the annihilation region, 
while the solid line corresponds to a simulation where the 
secondary nuclei are distributed homogeneously through- 
out the simulation volume. Note that photodisintegra- 
tion of 4 He was ignored in these calculations. We may 
thus compare the solid lines in panels (c) and (f) to the 
solid lines in panels (b) and (e). From this, it is evi- 
dent, that for annihilations occuring below T < 100 eV 
~ 80-90 % of all secondary nuclei are able to escape from 
the annihilation region and will thus survive subsequent 
annihilation. 



VI. DISCUSSION OF THE RESULTS 

A. Observational Constraints 

Any valid scenario for the evolution of the early Uni- 
verse has to reproduce the observationally inferred values 
of the light element abundances, which we will summa- 
rize in the following. The primordial 4 He mass fraction 
is commonly inferred from observations of old, chem- 
ically unevolved dwarf galaxies. Two distinct values 
are reported for the primordial 4 He mass fraction Y p — 
0.234±0.003 stat ±0.005 sys @ and Y p = 0.244 ±0.002 stat 
|Q . Very recently two of the pioneers of the field have 
determined the 4 He mass fraction in the Small Magellanic 
Cloud (SMC) by observing 13 areas of the brightest H II 
region in that galaxy: NGC 346 Hj. There observations 
are extrapolated to a value of Y p — 0.2345 ± 0.0026(lcr). 
While in excellent agreement with the above quoted lower 
value |42], it is in conflict with the higher value adver- 



tised in Ref. [Q and also not compatible with the cur- 
rently favored primordial 2 H determination. There are 
three claimed detections of ( 2 H/H) ratios at high red- 
shift from observations of quasar absorption line spec- 
tra. Similar to the case of 4 He, two conflicting values 
for the primordial ( 2 H/H) have been derived, 2 H/H = 



20 ±5 x 10" 



and 2 H/H = 3.39 ±.251 
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with stronger observational support for the low 2 H/H ra- 
tio. Using a new approach in analyzing the spectra, Lev- 
shakov, Kegel and coworkers reported a common value 
of 2 H/H = 4.4 ± .3 x 10~ 5 Q for all three absorption 
systems. Concerning the primordial 3 He abundance the 
situation is even less clear since only the chemically rel- 
atively evolved 3 He/H abundance in the pre-solar neb- 
ula is available, and the chemical evolution of 3 He is 
poorly understood. It is, however, reasonable to assume 
that the cosmic 3 He/ 2 H ratio is an increasing function 
of time. Whenever 3 He is destroyed in stars, the more 
fragile 2 H will certainly also be destroyed f4q| . Thus 
( 3 He/ 2 H) t > ( 3 He/ 2 H) p should hold for any time t after 

Big Bang nucleosynthesis. Given the pre-solar 3 He/ 2 H 
ratio |l9| we may impose the limit ( 3 He/ 2 H) p < 1. 

Only traces of 6 Li are produced in the framework of 
SBBN, and until very recently the observational data 
was very sparse. For this reason, 6 Li was not consid- 
ered to be a cosmological probe. With the confirmation 
of 6 Li detections in old halo-stars ]5(J and disk stars |5l]] 
on the level 6 Li/H ~ 7 x 10 12 this may change. Never- 
theless, since detection of such small 6 Li/H abundances 
requires operation close to the detection limits of current 
instruments, and possible stellar depletion of 6 Li is not 
well understood, the use of 6 Li as a cosmic probe may 
be controversial at present. In this light, we adopt a ten- 
tative upper bound on the primordial 6 Li abundance of 
about 6 Li/H ~ 7x 10~ 12 . This limit may be used to con- 
strain some non-standard BBN scenarios, which greatly 
overproduce 6 Li [jj6| . 

7 Li is inferred at a remarkably constant abundance of 
A( 7 Li)= log 10 ( 7 Li/H)- 12 = 2.238 ± 0.012 ±0.05 stat ||] 
in old POP II stars, referred to as the Spite-plateau. Nev- 
ertheless, stellar models which deplete 7 Li considerably 
have been proposed . Recently, it has been claimed 

that the primordial T7i abundance should even be lower 
than the plateau value 55 5(|. As in the case of 2 H, the 
7 Li abundance is however not useful to derive limits on 
our models. In the parameter range where the observa- 
tionally inferred 7 Li abundance yields are violated, limits 
derived from other elements are more stringent. 



B. Constraints on Antimatter Domains in the Early 
Universe 

In order to derive conservative limits on the amount 
of antimatter in the early Universe, we first discuss 
our results with respect to generally accepted observa- 
tional constraints. While there is currently a lively de- 
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bate about which of the two independent 4 He deter- 
minations reflects the primordial value, it seems rea- 
sonable to assume a 4 He mass fraction not lower than 
Y p 0.22 llH). No reliable limit on the 3 He abun- 
dance alone may be invoked, we therefore use the con- 
straint 3 He/ 2 H < 1 Q. These two values constitute 
our conservative data set as displayed in Fig. High 
antimatter fractions, Ra ^ 0.1 may only be consis- 
tent with the observationally inferred light element abun- 
dances if annihilation occurs close to weak freeze out, i.e. 
r A°° ^ 10 _1 cm. In this case, the weak interactions are 
still rapid enough to at least partially reproduce the an- 
nihilated neutrons and thus drive the n/p-ratio back to- 
wards the SBBN value. Antimatter fractions larger than 
Ra > a few 10 -2 on length scales r A ao > 10~ 4 cm result 
in an unacceptable low 4 He mass fraction, Y p < 0.22, 
which is indicated by the black shaded region in Fig. [n]. 
Even larger antimatter regions, r A 00 > 6 cm, annihilate 
at least partially via p He disruptions. Since the destruc- 
tion of only a minute fraction of 4 He leads to an obser- 
vationally unacceptable enhancement of the 3 He/ 2 H ra- 
tio, the limits on the allowed antimatter fraction in this 
regime may be as stringent as Ra ^ a few 10 -4 for length 
scales r A ao > 2 x 10 2 cm (dark grey shaded region in 
Fig. 0). Recently, Kurki-S uonio & Sihvola |ll|] derived 
similar limits based on the constraint 3 He/H < 10~ 45 . 
We feel that this choice is not optimal, given the large 
uncertainties in understanding the galacto-chemical evo- 
lution of 3 He. Furthermore, they find a production of 
3 He about a factor of 3 higher than in our study. 

If we employ the new and still slightly speculative 6 Li 
limit discussed above, we may significantly tighten the 
constraints on the amount of antimatter by requiring that 
pre-galactic production of 6 Li is not to exceed 6 Li/H < 
7 x 10 -12 . This leads to an improvement of the limit on 
Ra for late time annihilation, i.e. r A 00 > 6 cm, by up to 
two orders of magnitude. This is evident from the light 
grey shaded region in Fig. [Io| Nevertheless, due to the 
loophole of possible 6 Li depletion in stars, and due to the 
still preliminary nature of the 6 Li observations, this limit 
should be regarded as tentative at present. 

The limits derived from annihilations below T < 45 eV, 
corresponding to antimatter domain sizes of r A 00 > 
10 3 cm, have to be interpreted with care, since the pho- 
todisintegration yields in that regime are uncertain due 
to t he un known photon spectrum, as we discussed in 
Sec. [IV B| . But even if we completely ignore photodis- 
integration, meaningful limits due to direct production 
of 3 He (and subsequent 6 Li synthesis) via antiproton an- 
nihilation on 4 He may still be obtained. These limits 
are indicated by the dashed lines in Fig. M. Due to 
the increasing inefficiency of photodisintegration at low 
temperatures, both limits converge for large antimatter 
domain sizes. 

Note that the limits derived here on the basis of the 
3 He/ 2 H and 6 Li data should be similar in magnitude to 
limits on scenarios where antimatter is homogeneously in- 
jected into the plasma, for example by the decay of a relic 
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FIG. 10. Limits on the presence of antimatter in the early 
Universe. Parameter combinations within the black shaded 
region result in a 4 He mass fraction below 0.22, while in the 
dark grey shaded region the bound 3 He/ 2 H < 1 is violated. 
The excluded range is extended by the light grey shaded re- 
gion, if one adopts the tentative bound 6 Li/H < 7 x 10~ 12 . 
The dashed lines indicate the results when 4 He photodisin- 
tegration is ignored. Also shown are the usually weaker lim- 
its on the presence of antimatter from CMBR considerations 
(hatched region). 

particle after the nucleosynthesis epoch, since they rely 
on the production of secondary nuclei from 4 He disrup- 
tion and photodisintegration. Both processes are generic 
for scenarios with injection of antimatter. The competi- 
tion of annihilation within, and escape from, the annihi- 
lation region of the produced light isotopes is, however, 
particular to a scenario with individual domains. Escape 
of the annihilation products is inefficient for domain sizes 
between ~ 10 — 100 cm, corresponding to the "injection" 
of antimatter between temperatures ra 3 — 0.4 keV. In 
this regime more stringent constraints would apply to a 
homogeneous injection of the antimatter. Furthermore, 
the reduction of the n/p-ratio prior to 4 He synthesis and 
thus of the 4 He mass fraction also only applies to mod- 
els where antimatter is confined to well defined domains. 
Only in this case annihilation proceeds via baryon dif- 
fusion and thus the differential diffusion of charged and 
neutral baryons may provide an efficient sink for neu- 
trons. In contrast, a homogeneous injection of antimat- 
ter at temperatures above T4 He ra 80 keV (corresponding 
approximately to the scale r A 00 ra 6 cm in Fig. |l(J) may be 
constrained by an increase of Y p due to proton-neutron 
conversion induced by pion charge exchange ]32|| . 

For comparison, we have also shown in Fig. pi the lim- 
its on annihilation which may be derived from the upper 
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limits on distortions of the spectrum of the CMBR. The 
very precise CMBR data allows us to place constraints 
on the amount of non thermal energy input at redshifts 
below z«3x 10 6 . Each annihilation transforms about 
one half of the rest mass of the particles into electro- 
magnetically interacting particles, thus the limits given 
in Ref. ]58|] may directly be converted into a limit on 
i?A, which is indicated by the hatched region in Fig. [n]. 
Using the above conservative data set, we find stronger 
limits from BBN than the ones provided by the CMBR 
data for annihilations occuring at temperatures above 
?ann ^ 1 eV (7" A 00 < 10 5 cm), corresponding to a redshift 
of z > 4 x 10 3 . If we adopt the new 6 Li bound, the pres- 
ence of antimatter is more tightly constrained by BBN 
considerations, rather than by CMBR considerations, for 
the whole parameter range down to the recombination 
epoch at z « 10 3 . 



C. Upper Limit on Q,{, in Matter- Antimatter 
Cosmologies 

It is of interest to contemplate if a BBN scenario with 
matter-antimatter domains may reconcile the observa- 
tionally inferred element abundances with the theoret- 
ically predicted ones for a baryonic density exceeding 
the upper bound from SBBN, flbhf 00 < 0.02. Possi- 
ble alternative solutions to BBN which are in agreement 
with observationally inferred abundances for higher val- 
ues of SI5 have recently received renewed attention due 
to the results of the Boomerang and MAXIMA exper- 
iments |59| , |60|| on the anisotropies in the CMBR, which 
favor a baryonic density exceeding the SBBN value |3l| . 
In the standard BBN scenario, such a Universe suffers 
from overproduction of 4 He and 7 Li, and from severe un- 
derproduction of 2 H. In a scenario with annihilating an- 
timatter domains, there exist two possibilities to reduce 
the primordial He mass fraction to the observed value. 
Early annihilation, prior to 4 He synthesis, may reduce 
the n/p-ratio and thus the final 4 He mass fraction. Dur- 
ing late time annihilation 4 He nuclei may be destroyed 
via antiproton induced disruption and via photodisintc- 
gration. At first sight, the possibility to achieve obser- 
vationally acceptable 4 He mass fractions at high baryon- 
to-photon ratios looks promising. But upon closer in- 
spection, some severe shortcomings of such models arise. 
Scenarios at high net baryon-to-photon ratio and with 
annihilation prior to 4 He formation still overproduce 7 Li 
relative to the observational constraints. Furthermore, 
no additional source of 2 H exists in this model, which is 
thus ruled out. In the complementary case, where anni- 
hilation is delayed until after the epoch of 4 He synthesis, 
production of 2 H and 3 He due to disruption and photo- 
disintegration of 4 He results. Even though it is possible 
to find models where late time 2 H production may re- 
produce the observationally inferred value, the ratio of 
3 He/ 2 H will exceed unity. This is observationally un- 
acceptable. Further, such a scenario would produce 6 Li 



in abundance, which is most likely in conflict with re- 
cent observations. This remains true, even if we drop 
the assumption of a Universe in which the baryon- , or 
antibaryon-, to-photon ratio has initially the same value 
throughout the Universe, and furthermore allow the anti- 
matter fraction and domain length scale to take different 
values at different locations in space. Let us assume that 
the Universe consists of two different types of regions. In 
regions of type A with net baryon-to-photon ratio J7„ et , 
the antimatter fraction is high, i?A ^ 0.25, and mixing is 
effective between weak freeze out and T4 Ho . Irrespective 
of the exact value for r]^ ot , that region consists of protons 
only after the annihilation is complete. In region B, with 
net baryon-to-photon ratio f]^ ct , antimatter domains are 
larger, so they annihilate after the BBN epoch and light 
element synthesis may take place. If we further assume 
that region B is at high baryonic density, T]^ ei . 3> ^sbbn , 
the production of 2 H and 3 He is negligible prior to anni- 
hilation. Mass 2 and 3 elements will however be produced 
in the course of late time annihilation of p on 4 He. It is 
then easily feasible to find a ratio between the volumes of 
the two regions such that the average 4 He mass fraction 

0.25, 



Y p is diluted to the observed value of Y p 



(Yp) 



B 



1 



(35) 



where /a,b are the fractions of space occupied by the 
two different types of regions, respectively. The 4 He 
mass fraction converges to Y p « 0.36 for high baryonic 
densities, the required dilution factor is thus at most 
(Y p )b/Y p fn 1.5 in order to obtain Y p w 0.25. While 
the 4 He mass fraction may agree with the observational 
constraints for an arbitrarily large average baryon den- 
sity, we face the same problems with production of 2 H, 
3 He and 6 Li via late time annihilation as discussed above 
for a one zone model. Furthermore, in region B, 7 Li is 
produced well in excess of the observed values and the 
dilution by mixing with the proton-only zones of type A 
may not reduce the 7 Li abundance by more than a fac- 
tor of 1.5. We thus conclude that it seems difficult to 
relax the SBBN upper bound on Qj, by the existence of 
antimatter domains in the early Universe. 

A further result of our study is that the putative 
presence of antimatter in the early Universe may pro- 
vide some relieve for the tension between the lower of 
the two values for the primordial 4 He mass fraction, 
and the low 2 H determination, 2 H/H = 
In view of the recently reported value of 

, which co- 
incides with the low value, this discrepancy has received 
new attention. In a Universe at a comparatively high 
baryon-to-photon ratio of 77 w 5 x 10 -10 with an anti- 
matter fraction of a few 10 distributed on length scales 
smaller than 6 cm, the abundance yields for 4 He and 2 H 
may both be 'low' and thus the two observational con- 
straints mentioned above may be fulfilled simultaneously. 



Y p = 0.234 
3.39 xl0~ 5 _ 

Y p derived from observations of the SMC 
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VII. CONCLUSIONS 

In this work we have studied the mixing and subse- 
quent annihilation of antimatter domains in the early 
Universe during a period from a cosmological temper- 
ature of about 10 MeV, well above the epoch of weak 
freeze out, down to the formation of neutral hydrogen 
(recombination) at about 0.2 eV. Such distinct domains 
of antimatter may possibly arise in some electroweak 
baryogenesis scenarios (7|j|, as well as in other proposed 
solutions to the baryogenesis problem (for a review see 
e.g. |6|]). We have shown that the annihilation of anti- 
matter domains may have profound impact on the light 
element abundances. Depending on the time when anni- 
hilation occurs, we identify two main effects. Annihila- 
tion prior to the incorporation of all neutrons into 4 He 
results mainly in a reduction of the neutron-to-proton 
ratio, which determines the amount of 4 He synthesized. 
Such scenarios are thus constrained by the possible un- 
derproduction of 4 He. Even more stringent constraints 
on the antimatter-to-matter ratio may be derived if an- 
timatter resides in slightly larger domains and annihila- 
tion proceeds after the formation of 4 He. In this case, 
the dominant effect is the production of secondary ener- 
getic nuclei ( 2 H, 3 H, and 3 He), which may increase their 
respective abundances, but may also lead to the produc- 
tion of 6 Li nuclei. Further, energetic photons originating 
from the annihilation process may produce additional en- 
ergetic nuclei via the photodisintegration of 4 He. 

In a second aspect of our work we demonstrated that 
the presence of small amounts of antimatter, separated 
from matter within some length scale regime may, in fact, 
even improve the agreement between BBN theory and 
observations by reducing the amount of synthesized 4 He 
while leaving other light isotope yields basically unaf- 
fected. Finally, we argued that the SBBN upper bound 
on the cosmic baryon density fly, is very unlikely to be re- 
laxed in a scenario with annihilating antimatter domains 
in the early Universe. 
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APPENDIX A: 
STRUCTURE OF THE ANNIHILATION REGION 

In this appendix we will develop a detailed picture 
about the structure of the annihilation region at the 
boundaries between matter and antimatter domains. We 
will also illustrate why the numerical resolution of this 



thin layer is not essential for making fairly accurate pre- 
diction on the light-clement nucleosynthesis in an envi- 
ronment with matter-antimatter domains. 

In case of annihilation via neutrons, i.e. T ann > Tk Ho w 
80 keV, diffusion within the annihilation region is — ac- 
cording to the local baryon-to-photon ratio 77 ann — domi- 
nated either by magnetic moment scattering on electrons 
and positrons (ry ann < 10~ 8 ), or by nuclear scattering on 
protons (?7ann ^ 10 -8 ). In both cases, the typical scat- 
tering time for neutrons is much smaller than the anni- 
hilation time, 
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Here v b is a typical baryon thermal velocity, n™ 11 the 
antibaryon density in the annihilation region, ton the 
nucleon rest mass and the relevant cross sections are 



8 x 10 mb and a 



up 



2 x 10 mb (see e.g. 



Ref. [|20|]). Neutron scattering is thus always mo re prob- 
able than annihilation. Note that Eqs. (Al, A2) assume 
an electron density roughly equal to the photon density, 
n e ± w n-y, which is appropriate at early times i. e. before 
annihilation, when neutron diffusion is important. In 
the numerical computations, however, we follow the ex- 
act densities of the species. 

Annihilation via induced by protons diffusion occurs 
only in the keV era, where proton diffusion is limited 
by Thomson scattering of the electrons in the 'electron- 
proton system' off the ambient photons. Even though 
transport of the protons may now be controlled by hydro- 
dynamic expansion, the movement of the particles over 
the boundary and inside the annihilation region is still 
described by diffusion. Comparing the Thomson inter- 
action time with the annihilation time for protons (cf. 
Eq. ID, 
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we find that the scattering time scale is again much 
shorter than the annihilation time scale. In both cases, 
the width of the annihilation region l ann is thus given 
by the distance d(r ann ) nucleons can diffuse into the re- 
spective anti-region during their typical lifetime against 
annihilation r ann (cf. Eq. O), 
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We have included a factor of 2 to allow for diffusion of 
matter into the antimatter region as well as of antimat- 
ter into the matter region. The diffusion constant D can 
be taken to be constant over the lifetime against anni- 
hilation. In order to calculate r an n, we need to estimate 
the density in the annihilation region. We assume that a 
steady state between diffusion of baryon number into the 
annihilation region and annihilation of this baryon num- 
ber is established. The concept of a steady state is only 
appropriate for times somewhat shorter than the Hubble 
time, since the densities and diffusion constants vary with 
the expansion of the Universe. A typical baryonic density 
gradient some distance away from the annihilation region 
will however always be of the order of An& / (4(rHubbie) , 
with dfc(rHubbie) the diffusion length scale over one Hub- 
ble time. The difference in baryon density is given by 



An b = fi b 



1 , with fi b the baryon density far away 



from the annihilation region and n ann the baryon den- 
sity within the annihilation region. The baryon density 
in the annihilation region will typically be much smaller 
than fi b ; therefore we replace An b by fi b . This leads us 
to approximate the baryon number flux F b into the an- 
nihilation region by 



F b = D\7n b A ps D 
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The number of annihilations in a volume with surface A 
and width Z ann should then be equal to the flux of baryons 
into the volume, 
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As long as the diffusion length is considerably smaller 
than the size of the antimatter region, fib is equal to the 
initial matter density, hb = n net A°, where A is the ini- 
tial baryon overdensity and n net the initial average net 
baryon density (see Eq. |l|). We may now compute the 
baryon density in the annihilation region. Inserting the 



annihilation length Z ann , Eq. (A4), into Eq. (A6) and us- 



ing r ann = (cr ann u b n ann ) 1 yields 
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The baryon and antibaryon density within the annihila- 
tion region should be of the same magnitude, thus we 
finally obtain for the baryon density in the annihilation 
region 
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This may be written in terms of the local baryon over- 
density A ann in the annihilation region as 
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Interestingly, the overdensity in the annihilation region 
A ann is independent of the diffusion constant. We may 
now calculate the width of the annihilation region in our 
comoving units, 
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Using the relevant diffusion constants and annihilation 
cross sections and further assuming rj nct = 4 x 10~ 10 , we 
obtain 



Z^ n = 2 x 10~ 4 cm (A 



"1 1 ' ( -J—] 

MeV J 
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for annihilation via neutron diffusion and 



^0 = 1-3 x 10- 1 cm(A )- 1 / 3 



T 

keV 
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(All) 



(A12) 



for annihilation via proton diffusion. We have numeri- 



cally verified Eq. (All) for a scenario with antimatter 
regions of size rjf = 1.5 x 10~ 2 cm and an initial over- 
density A ps 10, corresponding to a matter-antimatter 
ratio of Ra = 0.9. In order to check the validity of the as- 
sumption of a steady state, we had to let the code evolve 
at least over the period of one Hubble time. The two 
snapshots of the neutron and antineutron overdensity, 
A„ and A fi (cf. Fig. [ll]) were obtained in a simulation 
which was started at T ss 20 MeV, and evolved down 
to T ps 10 MeV. The left panel shows the whole simu- 
lation volume, while the right panel is a zoom into the 
annihilation region of the same simulation. The resolu- 
tion is fine enough to describe (anti-) neutron diffusion 
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FIG. 11. Left panel: Snapshot of the neutron (full line) and the antineutron (dotted line) overdensity, A n and An at a 
temperature of T ~ lOMeV. This distribution was obtained with a high resolution simulation; the antimatter parameters were 
Ra = 0.9 and r A 00 = 1.5 x 10 -2 cm. Right panel: Zoom into the annihilation region. See text for discussion. 



within the annihilation region. We find an overdensity 
in the central region of A ann w a few 10~ 3 which may 
be compared to the above estimate, A ann « 3 x 10~ 3 . 
The width of the annihilation region is I^qq ~ 10 cm, 
following Eq. ( |All[ ). 

Since the two relevant processes — transport of par- 
ticles through their own region towards the annihilation 
region and diffusion within the anti-region — proceed on 
length scales which differ by orders of magnitude, it is 
very time-consuming to run simulations with the reso- 
lution necessary to adequately describe both processes. 
The numerical results presented in this work were ob- 
tained at a resolution which properly resolves the trans- 
port processes over the distance of order of the domain 
size, but does not resolve the diffusion within the anni- 
hilation region. This should however affect our results 
little, since the exact composition of the annihilation re- 
gion is not decisive for the final abundances. 

In case of annihilation before 4 He synthesis, the ex- 
act annihilation time is crucial for our results. Protons 
hardly play a role in case of early annihilation due to 
their very short diffusion length. The protons which are 
originally present in the annihilation region are quickly 
annihilated. Additional protons may not be transported 
into the annihilation region and their density profile re- 
mains frozen in. The annihilation region is thus pop- 
ulated by neutrons and antineutrons only, and further 
annihilation may only proceed via neutrons and antineu- 
trons. All particles which reach the annihilation region 
will inevitably be annihilated on a very short time scale 
compared to the transport time. Thus the time scale for 
annihilation of all antimatter is set by the transport of 
neutrons and antineutrons towards the annihilation re- 



gion, hence over considerably longer distances than the 
annihilation region, which are properly resolved. 

In case of annihilation after the disappearance of free 
neutrons at a temperature of T4 He w 80 keV, the domi- 
nant channels are pp and p He. The ratio of annihilations 
on either 4 He or on protons is important, since this ra- 
tio determines how many secondary nuclei, which arise 
in 4 He disruption, arc produced for a given antimatter 
fraction. This ratio depends again on the transport of 
the nuclei over the whole matter region into the anni- 
hilation region. The transport time scale may either be 
set by charged particle diffusion or by hydrodynamic ex- 
pansion. For both processes, resolution of the whole sim- 
ulation volume is important, but since again all nuclei 
which reach the annihilation region are inevitably anni- 
hilated, the spatial distribution of the nuclei within the 
annihilation region should be of negligible importance. 

The effect of not resolving the annihilation region is 
that matter and antimatter may travel further into the 
respective anti-region than is physically correct. But 
since in both cases discussed above the number of an- 
nihilations on a specific nucleus at a specific time is set 
by the transport processes, this lack of resolution should 
not be relevant. The relative independence of the results 
on the exact structure of the annihilation region is also 
evident by resolution studies given in App. [B|. 

Energetic secondary nuclei arising in the He disrup- 
tion process may only escape from the annihilation region 
and thus influence the final abundance yields if their stop- 
ping length is much larger than the annihilation region. 
The correct treatment of this effect is therefore indepen- 
dent of whether or not diffusion within the annihilation 
region is resolved. 
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APPENDIX B: NUMERICAL METHODS 

The task of performing detailed numerical BBN cal- 
culations was first fulfilled by Wagoner, Fowler & 
Hoyle fl|]. Jedamzik & Fuller [|o|Q developed an 
inhomogeneous BBN code to describe the evolution of 
subhorizon-scale baryon-to-photon fluctuations in the 
early Universe and the resultant modifications of the 
light element abundances. To this end, a Lagrangian 
grid of zones was introduced in which the various ther- 
modynamic quantities and the nuclear densities may de- 
viate from the respective horizon average value. The 
BBN network is coupled to all relevant hydrodynamic 
processes, such as diffusion of baryons, photon diffusive 
heat transport, neutrino heat transport and late-time hy- 
drodynamic expansion of high-density regions. The nu- 
clear reaction network and the thermodynamic evolution 
of the homogeneous radiation background is treated as 
in an updated version of the original BBN code |^3|,|l9| . 
Baryon diffusion and incorporation of baryons into nu- 
clei proceeds on fast time scales, it is thus necessary to 
treat baryon diffusion implicitly. Further, neutrino and 
photon heat transport and hydrodynamic processes are 
included in the code. 

We extended the inhomogeneous code |2(],|d]] to in- 
clude antielements and adapted it to the present prob- 
lem. We use a set of concentric spherical shells to de- 
scribe the distribution of matter and antimatter. The 
number of zones in the simulations has to be chosen such, 
that the spatial resolution of the volume is sufficient to 
adequately describe the relevant physical processes. It 
turned out that 30 zones are sufficient to resolve diffu- 
sion of the nuclei and obtain reasonable convergence in 
the final abundances (see Fig. [l2]). The region where 
most of the annihilations occur, on the other hand, has 
not been resolved in our simulations. The results should 
however be fairly insensitive to the actual structure of 
the thin annihilation layer, as we discuss in App. [A|. 

The whole procedure of solving the nuclear reaction 
matrix and the treatment of the hydrodynamical pro- 
cesses is included in a second order Runge-Kutta scheme, 
i.e. it has to be done twice per time step. The results 
of the two Runge-Kutta steps are averaged at the end of 
each time step. Since not only the densities within each 
zone, but also the radii of the zones may change, this has 
to be done very carefully in order to minimize errors. An 
adequate independent check of the numerical simulations 
is the achieved accuracy in baryon number conservation. 
Generally, the baryon number is conserved on the level 
of ANb/Nb < C(10 -6 ) for those regions of the parame- 
ter space relevant for the derivation of our limits. Only 
for large antimatter fractions, Ra > 0.5 on length scales 
r A°° ~ 10 cm, numerical conservation of baryon number 
deteriorates. It seems however very unlikely to find an 
observationally acceptable scenario for combinations of 
the parameters in that range such that the simulation of 
such scenarios is of little interest. 
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FIG. 12. Light element abundances obtained in simulations 
with three different combinations of the antimatter parame- 
ters (see legend). The number of zones was varied between 
20 and 70 to check for convergence of the results. 



The escape of 6 Li and the mass three nuclei from the 
annihilation region was treated as follows. We keep track 
of the number of 3 H, 3 He and 6 Li nuclei produced in 
p 4 He annihilations during a single time step. Using the 
stopping length for these nuclei according to Eq. (Eq) we 
may calculate the fraction of those thermalized within 
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the matter region by geometrical considerations. This 
fraction is then added to the number density of the re- 
spective nuclei. 
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